Emission of photon pairs at discontinuities of nonlinearity in spontaneous parametric 

down-conversion 
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In order to fulfil the continuity requirements for electric- and magnetic-field amplitudes at discon- 
tinuities of x' 2 ' nonlinearity additional photon pairs have to be emitted in the area of discontinuity. 
Generalized two-photon spectral amplitudes can be used to describe properties of photon pairs gen- 
erated in this process that we call surface spontaneous parametric down-conversion. The spectral 
structure of such photon pairs is similar to that derived for photon pairs generated in the volume. 
Surface and volume contributions to spontaneous down-conversion can be comparable as an example 
of nonlinear layered structures shows. 



I. INTRODUCTION 

The generation of second-harmonic field at a bound- 
ary between two homogeneous media with different val- 
ues of x nonlinearity has been addressed for the first 
time more than 30 years ago [3, [§] • This weak effect has 
been discovered when second-harmonic generation with 
considerable phase mismatch has been investigated. The 
surface second-harmonic field occurs here naturally and 
assures the fulfilment of continuity requirements for the 
tangential components of electric- and magnetic-field vec- 
tor amplitudes that stem from Maxwell's equations. In 
more detail, a fundamental field creates a nonlinear po- 
larization at second-harmonic frequency at the nonlinear 
side of the boundary. This polarization generates two 
surface second-harmonic fields, one in forward direction, 
one in backward direction. As a consequence two differ- 
ent second-harmonic fields propagate inside a nonlinear 
crystal. They differ in their wave vectors. The first (and 
usual) field originates from the volume nonlinear polar- 
ization and its local wave vector is twice the wave vector 
of the fundamental field. On the other hand the wave 
vector given by index of refraction at second-harmonic 
frequency characterizes the second-harmonic field arising 
at the boundary and propagating freely through the crys- 
tal. Experimental evidence of these effects can be found, 
e.g., in 0. Surface second- harmonic generation pumped 
by ultrashort pulses has been analyzed in Q. Deep un- 
derstanding of this effect can be reached when studying 
this process in a nonlinear medium with negative index 
of refraction [4J. In Ref. 0], completely numeric ap- 
proach based on the solution of nonlinear Maxwell's equa- 
tions has been adopted contrary to the original analytical 
and approximate approach in [l| demonstrating the rich- 
ness of physical effects included implicitly in Maxwell's 
equations. We note that also inhibition of absorption in 
highly phase-mismatched volume second-harmonic gen- 
eration has been observed 

The above described effects are valid for nonlinear 
parametric (three-mode) interactions in which a large 
number of material states far from resonance participate. 
On the other hand resonant second-harmonic generation 
meadiated by resonant surface states has been widely 



studied for many materials (see, e.g., in p, Q) and has 
become a useful tool for surface diagnostics at present. 
We note that also entangled photon pairs generated in 
parametric down-conversion can be converted resonantly 
into plasmons at material surfaces even with the preser- 
vation of polarization entanglement @. 

The question arises whether nonresonant surface ef- 
fects can occur also in the quantum process of spon- 
taneous parametric down-conversion (SPDC) 0. In 
the case of second-harmonic generation, the presence of 
macroscopic classical nonlinear polarization at frequency 
2w is crucial. On the other hand, there is no macroscopic 
classical nonlinear polarization at the signal- or idler-field 
frequencies. However, quantum nonlinear polarization 
occurs in SPDC at these frequencies and is responsible 
for photon-pair generation at a boundary. 

SPDC is described by an appropriate momentum oper- 
ator that is constructed in the framework of energy-flux 
quantization [Tl|, [H, EH ■ We note that the inclusion of 
all fields occurring during the propagation, i.e. forward- 
as well as backward-propagating fields, is necessary to 
keep consistency of the approach. 

The article is divided as follows. A model of surface 
SPDC in case of a homogeneous nonlinear crystal is de- 
veloped in Sec. II. Determination of physical quantities 
characterizing photon pairs is described in Sec. Ill using 
generalized two-photon amplitudes. Sec. IV gives a gen- 
eralization to nonlinear layered structures. Conclusions 
are drawn in Sec. V. 



II. MOMENTUM OPERATOR AND FIELDS' 
CONTINUITY AT THE BOUNDARIES OF A 
NONLINEAR CRYSTAL 



In this section, we first pay attention to the volume 
nonlinear interaction, then study the problem at the in- 
put and later at the output boundaries and finally add 
the obtained expressions describing the photon-pair gen- 
eration. We note that a simplified model has been pre- 
sented in 1141 . 
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Volume interaction 



The following interaction momentum operator Gi nt is 
appropriate for the process of SPDC [t| [ll| : 



Gi n t(z) = <ie A dt ^2 d-i,ap 

a,P,<y=F,B 



X 



4+ ) (z,^t)(z,t) J B} (3 '(z,<) + h.c. , (1) 



where Ep^ are positive-frequency parts of the (linearly 
polarized) pump-field electric-field amplitudes whereas 

Eg a {e\ ' ) stand for negative-frequency parts of the 
signal- (idler-) field electric-field amplitude operators. 
Subscript F (B) refers to a field propagating forward 
(backward), i.e. along the +z (—z) axis. Symbol eo 
means permittivity of vacuum, d 7iQ ^ are effective non- 
linear coefficients, A is transverse area of the fields, and 
h.c. replaces the hermitian-conjugated terms. We have 
assumed a scalar model for the interacting fields for sim- 
plicity. However, a generalization to the vectorial model 
is straightforward because of the applied first-order per- 
turbation approximation. Using a spectral decomposi- 
tion of the interacting fields, 

E ml ( z > t ) = "7S= / duj ™ E ml ( z ' exp(-MJ m *) , 

V 2lT J 

m = p,s,i, a — F,B, (2) 

the momentum operator G; n t in Eq. (fT]) can be recast 
into the form 



Gint (?) 



x E 



27T 



dujp I dio s 



duji 8(u) p — lo s — u>i) 



X 



[Eg (z, u p )E£(z, lo s )e\; ] (z, Ui ) + 



h.c. 



(3) 



where the 5-function expresses conservation of energy 
for monochromatic waves. The signal and idler spectral 
electric-field amplitude operators Em] can be expressed 
in terms of creation operators a ma introduced such that 
a m o mct gives the photon-number density in mode m a at 
a given frequency: 



huj r/ 



aL(z,u m ), 



2e cAn m (uj m ) 
m = s,i, a — F, B 



(4) 



Symbol n m stands for an index of refraction of field m. 

Spatial evolution of optical fields is determined by the 
solution of Heisenberg equations [ll| [l6| for field opera- 
tors denoted as X: 



dX{z) 
dz 



G(z),X(z) 



(5) 



G(z) 
G (z) 



Go(z) + G int (z), 

h 



E E 

m—s,i a—F.B 



(6) 



the interaction momentum operator Gmt is given in 
Eq. ([3]). The momentum operator Go introduced in 
Eq. ^ describes free-field evolution. Symbol k ma stands 
for a wave vector of mode m a at frequency u> m ; k mp = 
k m , k mB = —k m , k m > 0. Symbol h denotes the reduced 
Planck constant. 

In order to determine the electric-field amplitude oper- 
ators at the output of the nonlinear crystal, the Heisen- 
berg equations for operators d mQ (z, w m ) (m = s,i, a = 
F, B) have to be solved: 



da Sa (z,uj s ) 
dz 

+ E 

x exp[ik p ^(u} s + uJi)z}dl (z,0Ji), a 



ik Sa (uj s )a Sa (z,uj s ) 
dw i g lta f 3 (u) s ,Lo l )E { +X(),u) s +u)i) 

F. B. (7) 



Equations ([7]) have been derived assuming equal-space 
commutation relations [HJ. We note that fields propa- 
gating along the —z axis have negative wave- vectors in 
the definition of the free-field momentum operator Go 
in Eq. ([6]). A more detailed and rigorous formulation of 
the dynamics of counter-propagating fields justifying this 
approach can be found in [15(. The coupling constants 
57, a/3 occurring in Eq. ([7]) are given along the expression: 



2id 



27T7l s (uJ s )ni(uJi) 



(8) 



Equations for operators di F and &i B can be derived from 
Eq. by the formal substitution s <-> i. 

Solution of Eq. ([7]) valid up to the first power of g can 
be obtained in the form: 

a Sa (z,Lu s ) = e-xp[ik Sa (uj s )z] [a Sa (0, u) s ) 



E 

0,y=F,B 



duiB 7 , a p(z, aj a ,Ui)al (0, w 4 ) 



a 



F,B, 



(9) 



where 



x exp[iAkj^ a p(uj s ,uJi)z/2] 
xzsmc[Ahy t a/3(w a ,u/i)z/2]-, (10) 
Afc 7iQ/3 (w s , Wi) = kp^(bj a +uji)-k aa (uj a )-ki g (ui); 

a,0,>r = F,B; (11) 

sine (a;) = sin(x)/a;. We note that the restriction to 
the first power of g in the formula in Eq. ^ has jus- 
tified the use of the approximate formula a mQ (z,aj m ) = 
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exp[ik rrla (uj m )z]a ma (0, ui m ) below the integral over the 
frequency Wj. The solution in Eq. © describes SPDC 
originating in the volume of nonlinear crystal. 

The solution obtained for annihilation operators 
a Sa (z, u) s ) as written in Eq. ([9]) provides the following 
expressions for the positive-frequency parts of electric- 



(E^ ) and magnetic-field (H^') amplitude operators 



2e cAn s {u) s 



■ exp[ik Sa (w s )z] 



,(o,w.)+ ]T 

P,-y=F,B 



dhJiB lta p{z, Lo s ,uji)a\ (0, u>i) 



(12) 



(13) 
(14) 



he 



2/j, uj s A 



n s (uj s ) 



x ^2 / du>ig Jta p(ij a ,u)i 

x exp[ifcp T (w s + Wi)z] exp[-ife ip (a; i )2]at p (0,a;i) ) 

a — F,B. (15) 

Equations ([1151 — fT5"|) for the magnetic-field amplitude op- 
erators H Sa (a = F 1 , -B) have been derived assuming po- 
larization of electric-field amplitudes E Sa along the +x 
axis and, consequently, polarization of magnetic-field am- 
plitudes H Sa along the +y axis. The Maxwell equa- 
tions then provide the following formula H^\z,lu s ) — 

—i/(u> s iJ,o)dEs^ (z,u) s )/dz, where fio stands for perme- 
ability of vacuum. The magnetic-field amplitude op- 
erators Ht£\z,u) a ) have been decomposed in Eq. (fT3")) 
into two parts; the amplitude operators Hs^ (z, uj s ) 
are linearly proportional to the electric-field amplitude 
operators E^ (z,ui s ) whereas the amplitude operators 
,uj s ) are of purely nonlinear origin. The am- 
plitude operators r are not taken into account in 
the usual derivation of Fresnel's relations that assumes 
linear media. Correct inclusion of these amplitude op- 
erators into the continuity considerations at a boundary 
results in additional contributions to the nonlinear pro- 
cess. 



B. Input boundary 

Let us take deeper attention to the problem of con- 
tinuity of electric- and magnetic-field amplitudes at the 



boundaries of the nonlinear medium. First we pay at- 
tention to the input boundary (z — 0) and consider 
the signal field. Four electric- and magnetic-field am- 
plitudes are involved in the continuity requirement: am- 
plitudes Es£(0) and H^(0) of the forward-propagating 
field impinging on the boundary from the outside of non- 
linear medium, amplitudes Eg J (0) and HjQ (0) leaving 
the boundary outside the nonlinear crystal, amplitudes 
E Sg (0) and H SB (0) of the field coming to the bound- 
ary from the nonlinear crystal, and amplitudes E SF (0) 
and H SF (0) leaving the boundary and propagating in- 
side the nonlinear crystal [see Fig. {I}]. Because the 
magnetic-field amplitudes H Sp (0) and if (0) defined in 
the nonlinear crystal have also nonlinear contributions 
if s n J r (0) and Hfj^iO) described in Eq. ([15]) additional 

(surface) amplitude corrections SE 8p (0) and 5Ei° B (0) to- 
gether with SH SF (0) and 5 Hag (0) are needed to fulfil the 
continuity requirement. The surface amplitude correc- 
tions naturally occur in the fields that leave the bound- 
ary which is a consequence of spatio-temporal consider- 
ations that are suppressed to certain extent in our one- 
dimensional model. 

The requirement of continuity of projections of electric- 
and magnetic-field amplitudes to the plane of the input 
boundary leads to the following equations: 

Ei°J (0) + E<® (0) + SE(°] (0) = E SF (0) + 8E SF (0) 

+E SB (0), (16) 
(0) + (0) + (0) = Hf F (0) + if » F Fr (0) 

+SH SF (0) + H^ B (0)+H^(0). (17) 

Derivation of the usual Fresnel's relations (in linear me- 
dia) [22j is based on the fulfilment of the following equa- 
tions: 

£(°)(0) + £j°)(0) = E SF (0) + E SB (0), (18) 
Hi°J(0) + Hi o B \0) = <(0) + <(0). (19) 

Comparison of Eqs. (fT6"|) and (fT5|) together with Eqs. (fTTj) 
and (fT9"|) results in two algebraic equations for the surface 
amplitude corrections of the fields leaving the boundary: 



5E2»(p) = 5E SF (0), 



6H SF (0) + H^(0). 



(20) 
(21) 



Alternatively and more conveniently, the amplitude 



corrections 6Ei° B and 6Hl B ' of the field outside the non- 
linear crystal can be formally included into the equations 
giving Fresnel's relations. This can be done if we intro- 
duce fictitious amplitude corrections SE SB and SH SB of 
the field impinging on the boundary from its nonlinear 
side. Such corrections give, after transformation at the 
boundary using Fresnel's relations, the required ampli- 
tude corrections 5E S ° B and SH^J [see Fig. [TJ . Then we 
have: 

E S °J (0) + \Ei°J (0) + 5E<® (0)1 = E SF (0) 



r(0) 
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B(°)(0) + iS£'°l(0) 



E.JO) +SE.JO) 



E, B (L) + 6E, B (L) EW(L) 



E„(0)+SE„(Q) 



E„(L) +SE„(L) Eg)(L) + SEg)(L) 



FIG. 1: Scheme showing electric-field amplitudes E and their 
surface corrections SE at the input (z = 0) and output 
(z = L) boundaries of a nonlinear crystal. Superscript (0) 
[(1)] denotes amplitudes in front [beyond] the nonlinear crys- 
tal. Amplitude corrections SE SB (0) and SE 3F (L) written 
in frame-boxes do not exist in the real nonlinear medium; 
they replace the effect of real amplitude corrections 8E 3F (0), 
SE^(0), 8E%(L), and SE SB (L). 



+ [E SB (0)+6E SB (0)} ) (22) 

+ [HJl(0)+SH SB (0)]. (23) 

Equations p^|) and (jTTJ) are then fulfilled provided that 
the following two algebraic equations for surface ampli- 
tude corrections of the fields inside the nonlinear crystal 
are valid: 

= 6E SF (0) - 6E SB (0), (24) 
= H*? (0) + 6H SF (0) + - SH tB (0). 

(25) 

The positive-frequency parts of surface amplitude- 
correction operators 8E ma and SH ma are defined simi- 
larly as the corresponding amplitude operators E ma and 
H ma in Eqs. and (|14p using operator corrections 5d ma 
to the annihilation operators a ma : 



6E^(z,u) m ) = 



huj r/ 



2e cAn m (uj m ) 



5a ma (z, u m ),(26) 



m = s,i, a = F, B. 



(27) 



Substitution of Eqs. |T5]), ((26j) , and d2TJ) into 
Eqs. (f24|) and ([25| gives two algebraic equations for 
the annihilation-operator corrections 5d SF (0,uj s ) and 
Sa SB (0 1 uj s ): 



Sa SF (0,u> s ) - 5a aB (0,oj a ) = 0, 
ik SF (uj s )Sa SF (0, uj s ) - ik SB (u s )Sa SB (0, uj s ) 

a,p,~ i =F,B' ! 

xswcw. + woaj (o,wi) = 0. 



(28) 



(29) 



Solution of Eqs (|28|) and (j29|) finally gives the expres- 
sions for annihilation-operator corrections 5a SF and 5a SB 
at the input boundary: 



Sa SF (0, u> s ) 



2k s {ui s ) 



^2 / dcJ i9~/,a(}(u s ,UJi) 



a,/3,j=F,B ' 

<5a SB (0,w s ) = Sa SF (0,uj s ). 

C. Output boundary 



(30) 
(31) 



Fields at the output boundary of the nonlinear crystal 
can be analyzed similarly as at the input boundary. Here, 
the continuity of projections of electric- and magnetic- 
field amplitudes at z = L (L denotes the crystal length) 
to the plane of the boundary gives two equations: 



E SP (L) + E SB (L) + SE SB (L) 

= e£(L) + 6EW(L)+eU(L), 



(32) 



HZ (L) + Hlf(L) + HZ (L) + HfZ(L) + SH SB (L) 

(33) 



hU(L)+SHW(L)+hW(L). 



Amplitudes E^J (L) and Hs F (L) describe the field out- 
side the nonlinear crystal whereas amplitudes Es b \l) 
and Hg B {L) refer to the field impinging on the output 
boundary from its linear side (see Fig. |T|) . The amplitude 
corrections SE^J (L) and SHi F \L) can be formally in- 
cluded into the equations that express Fresnel's relations 
provided that fictitious amplitude corrections 5E SF (L) 
and SH SF (L) are introduced. Motivation for this step is 
the same as in the case of input boundary: we want to 
have corrections only inside the nonlinear crystal. Wc 
then have: 

[E SF (L) + 6E SF (L)} + E SB {L) = 

'eW(L) + 6EW(L)} +e£(L), (34) 

[hJ;(l) + 6h Sf (l)]+hZ(l) = 

'hW(L) + 6HW(L)]+hW(L). (35) 



r(i), 



Comparison of Eqs. ([M]) and (03 with Eqs. (J32I) and 
(1331) results in two algebraic equations for the surface 
amplitude corrections inside the nonlinear medium: 



-SE SF (L)+5E SB (L) = 0, 
H«Z (L) - 6H SF (L) + Hf^(L) + 6H SB (L) = 0. 



(36) 



(37) 



Two algebraic equations for the annihilation-operator 
corrections Sa SF (L,Lu s ) and Sa SB (L,uj s ) can be derived 
from Eqs. (|3"6"|) and (f3"7) using the expressions in Eqs. 
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JT5L (ESI), and (HTJ): 



Sa SF (L, u s ) + 5a SB (L,u s ) = 0, 



(38) 

-ik sp (u s )Sa SF (L, u s ) + ik SB (u s )Sa SB (L, u s ) 

a,/3,j=F,B 

x E^> (u s + Ui) exp[iA: P7 (u s + u t )L] 

x exp[-ifcj (3 (a; i )i]aj 9 (0,a; < ) = 0. 

(39) 

Equations (|38|) and (139)) can be solved leaving us 
the expressions for annihilation-operator corrections 
5a SF (L,u s ) and Sa SB (L, u s ) at the output boundary: 



Sa SF (L,u s ) = 



2k s (u> s ) 



xE^'(u s + Ui) ex.p[ik Py (u s + u t )L] 
x exp[-ik i0 (Ui)L]al (O,Ui), (40) 
5a sg (L,u s ) = Sa SF (L,u s ). (41) 

D. The whole nonlinear crystal 

The overall solution for operators a Sa {L,u s ) valid up 
to the first power of g has three nonlinear contributions: 
the first comes from the input boundary, the second from 
the volume and the third from the output boundary. The 
overall solution can be written as follows: 

a Sa (L,u s ) = a^ c (L,u s ) + £ 

d^i ^,af3( L , Us, u t )a™ cf (L, Ui), a = F,B. (42) 

Operators a l ™°(L, u m ) (m = s,i, a — F,B) oc- 
curring in Eqs. (|42|) and (|43|) below and expressed 
at the crystal end correspond to free-field evolution 
(i.e., without photon-pair generation) inside the crys- 
tal. They are determined by the formula a?™°(L,u m ) — 
exp[ifc„ lci (u m )L]a^(0, u m ). Functions T^ a p are defined 
in Eqs. (|4H- t4"Hjl bellow. 

Now we pay attention to the idler fields and use sym- 
metry between the signal and idler fields. The idler-field 
electric- (magnetic-) field amplitudes are assumed to be 
polarized along the +y (—x) axis. The requirement of 
continuity of electric- and magnetic-field amplitudes at 
the input and output boundaries leads to the solution 
for idler-field operators ^(L,^) in the form: 

a i/3 (L,Ui) = a l * cc (L,u l ) + ^ 

a,~/=F,B 

[ du s T\^{L, u s , Ui)a^(L, u s ), (3 = F, B. (43) 



The functions T s and T % introduced in Eqs. (|42f and 
(|43l) can be decomposed into volume (jF vo1 ) and surface 
(J s ' surf , jr^urf) contributions: 



= t: 



vol 



m,surf 



7, a/3 — J 7, a/3 T" ^'j.aP ' (44) 

F^lpiL^s,^) = g ltO , (u s ,Ui)E^(u s +Ui) 
x exp[ifcp T (w s +Ui)L] exp[-iAfc 7!Q/3 (a; s , Wi)L/2] 
xLsinc[A/c 7jQ/ 3(w s , Ui)L/2], (45) 

FT£p Tt {L,u s ,Ui) = — ■ — rg yta p(u} s ,Ui)E^(u B + u{) 

Km \U m ) 

x {exp[ifc Sct (w s )£] eiq>[iki p (uJi)L} 
- exp[ik Py (u s + Ui)L]} , 

m = s,i; a,/3,j — F,B. (46) 

In deriving Eq. (|46[) we have assumed g 7 ,F/3 = 9j,b/3 and 
g-y, aF = g~i,aB- The functions F m > sulf describing surface 
contributions disappear in the limit L — > 0, i.e. the sur- 
face contribution from the output boundary completely 
compensates that from the input boundary. Comparison 
of the expressions in Eqs. (|4"5| and (|4*6")l reveals a simple 
relation between the volume and surface contributions: 

•7-m,surf / t- \ . . f , 

• / " 7 ,q/3 \L>,U s ,Uj) Ak ltaj3 {UJ s ,UJi) 



^(L,u s ,Ui) 



v 7 m Q/3 K,^), 



x,p, 1 = F,B. (47) 



The structure of surface contribution is formed by mutual 
interference of fields generated at the input and output 
boundaries. At a boundary, only the energy conserva- 
tion restricts properties of an emitted photon pair. Such 
photon pair thus has a rich internal spectral structure be- 
cause phase-matching conditions do not apply. It is then 
the mutual interference of fields coming from the input 
and output boundaries that gives conditions similar to 
those of phase matching naturally found in the volume 
interaction. 

In standard bulk sources of photon pairs that are typ- 
ically several mm long, the interaction among forward- 
propagating pump, signal, and idler fields is important. 
In this case, the commonly-used formalism for the de- 
scription of SPDC [using a two-photon spectral ampli- 
tude &(u s ,Ui)] can be applied (see, e.g., [13, El, fl3|)- 
The surface contributions can be involved in this formal- 
ism if the following formal substitution is done: 



<S>(o> s ,cJi) 



1 + V 



F.FF 



{u s ,Ui) 



x yfi + V l FiFF {u s ,u t ) <5> m \u s , Ui), (48) 



where the two-photon amplitude & vol (u s ,Ui) charac- 
terizes the usual volume contribution. Under the 
usual condition of phase-matched nonlinear interaction 
[Akppp(u®,u®) = 0] the surface contributions occur 
only at spectral tails and are negligible. 
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Surface contributions have typically broader spectra 
compared to the volume contributions. If these spec- 
tra are not filtered in an experimental setup they lead 
to sharper features of two-photon temporal amplitudes 
JP s < surf (T s ,T 4 ) and ^i, S uxf( Ts)T .) [ for their definition, see 

Eq. (f56|) below]. As documented in Fig. [2ja, for a 
BBO crystal, the surface two-photon temporal ampli- 
tudes jF m ' surf (m = s, i) attain large values in the vicin- 
ity of the input and output boundaries. On the other 
hand the two-photon temporal amplitude JF vo1 of the vol- 
ume contribution has roughly the same values along the 
whole nonlinear crystal. If the volume SPDC is strongly 
phase mismatched its two-photon temporal amplitude 
.F™ 1 resembles that of the surface SPDC (see Fig. Eb). 
This means that only values of the two-photon temporal 
amplitude ^ vo1 characterizing photon pairs born in the 
vicinity of crystal edges are higher. Destructive interfer- 
ence inside the nonlinear crystal prevails in this case and 
suppresses photon-pair emission. Because the signal- and 
idler-field spectra are very broad for phase-mismatched 
interaction the shapes of two-photon temporal ampli- 
tudes T for the volume and surface interactions are sim- 
ilar. We also note that the temporal widths of peaks of 
amplitudes F s (t s ) in Fig. [2] are broader in the area that 
corresponds to the beginning of the crystal compared to 
those coming from the crystal end (and occurring around 
t s k s) because of intermodal dispersion faced by a 
photon-pair as it propagates through the crystal. 

The fact that only photon pairs around the bound- 
aries are generated in the strongly phase-mismatched in- 
teraction resembles the behavior of the second-harmonic 
field in the process of strongly phase-mismatched second- 
harmonic generation [5j. Here, the pulsed second- 
harmonic field propagates below the fundamental pulsed 
field (they have the same group velocities) along the crys- 
tal and does not feel any absorption |5|. This can be in- 
terpreted so that the second-harmonic field at the crys- 
tal output is generated only in the vicinity of the output 
boundary. 

Two-photon temporal amplitudes are not experimen- 
tally accessible but certain information about their shape 
can be reached when measuring coincidence-count 
interference rates R n [defined in Eq. (|58|) below] in a 
Hong-Ou-Mandel interferometer. Whereas a triangular 
dip is typical for phase-matched volume SPDC, two side 
dips with reduced visibility (around 0.5) and one central 
peak occur in the coincidence-count rate R n for surface 
SPDC as a consequence of the shape of two-photon tem- 
poral amplitude with two peaks (see Fig. [3]). We note 
that the profile of coincidence-count rate R n of strongly 
phase-mismatched volume SPDC is similar to that appro- 
priate for surface SPDC. Alternatively, sum-frequency 
generation of the signal and idler fields can be used to 
experimentally 'scan' the shape of two-photon temporal 
amplitude [2(J, Hl[ . The intensity field arising from sum- 
frequency generation should have two peaks depending 
on the mutual temporal delay of the signal and idler pho- 
tons in the studied phase-mismatched case. These peaks 



b) 




-30 30 60 90 120 150 180 210 240 

t s (l(T 14 s) 

FIG. 2: Cross-section of the absolute value of two-photon 
temporal amplitude F a (T a ) = \T% FF (r 3 , 0)| for the surface 
(b = s, surf, solid line with *) and volume (b = vol, solid line) 
contributions assuming cw pumping. The values of l-T^'^^l 

are 5000 (10) times magnified with respect to these of 
in case a (b). The curves are appropriate for a BBO crys- 
tal 5 mm long and collinear type-II interaction at the pump 
wavelength of 400 nm and signal and idler wavelengths 800 nm 
(frequency filters 30 nm wide (FWHM) are used). The crystal 
optical axis declines by 42.35 deg (perfect phase matching) (a) 
and 80 deg (b) with respect to the axis of fields' propagation. 



correspond to two boundaries. 

Contrary to the case of usual nonlinear crystals, surface 
SPDC cannot be neglected in nonlinear layered media 
(see Sec. IV later). 



III. QUANTITIES CHARACTERIZING THE 
EMITTED PHOTON PAIRS 

The inclusion of surface contributions to the expres- 
sions for operators a SF , a SB , cii F , and &i B requires a 
generalization of the usually used formulas for the de- 
termination of physical quantities characterizing photon 
pairs. For simplicity we pay attention to photon pairs 
with both photons propagating forward and define op- 
erators a m (u m ) ee t m (u m )a mp (L,uj m ) and Em\r m ) = 
t m {^m)Eml{L,T m ), where the coefficients t m describe 
the amplitude transmissivities at the output boundary 
(m = s, i). We have t m = 2n m / [n m + rim) according to 

Fresnel's formulas, n m (rim^) means an index of refraction 
of (beyond) the nonlinear medium (m — s,i). 
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1.5 



1.0 



0.5 



0.0 




pairs: 



r: (10- 14 s) 



FIG. 3: Normalized coincidence-count rate R n in Hong-Ou- 
Mandel interferometer as a function of relative time delay t; 
is shown assuming volume (solid line) and surface (solid line 
with *) SPDC. Values of parameters appropriate for Fig. 2a 
are used. 



The joint signal-idler photon-number density n(u) s , cjj) 
giving the number of emitted photon pairs with a signal 
photon in the unit interval around frequency lo s and its 
idler twin in the unit interval around frequency u>i at the 
medium output plane is defined as follows: 

n{uj s ,LUi) = ^ al(u s )a s (u s )al(u}i)ai(u}i) +h.c. ^ /2, 

(49) 

Symbol ( ) means averaging over the initial signal- and 
idler-field vacuum states. Substitution of Eqs. (|42[) and 
(H2J) into the definition of the joint signal-idler photon- 
number density n in Eq. (|49[) results in the formula: 



n(cj s ,u)i) = Re{F s *(ui s ,uj l )J :t (u) s ,uji)}. 



(50) 



The functions J- s and T % in Eq. (|50[) include transmission 
through the output boundary: 

T m {u: s ,uji) =t s (uj s )t i (uj i )T™ FF (u) s ,uji), m = s,i. (51) 

In Eq. ([S1I)) . symbol Re denotes the real part of an argu- 
ment. We note that expressions containing factor S 2 (uj p ) 
occur in formulas like that written in Eq. (|50[) for cw 
pumping. This factor has to be replaced using the for- 
mula S 2 (ojp) — IT I (2-K)b(io v ) in this case where 2T means 
the detection-interval length 

Intensity spectrum 5 s (w s ) of, e.g., the signal field can 
be easily derived from the joint photon-number density 



S s (uj s ) = huj s / dwjn(w 8 ,Wi) 
Jo 

/>oo 

= huj s / cL'jRe{^' s *(w s ,w i )^ i (a; s ,w i )}. 
Jo 

(52) 

Photon flux of, e.g., the signal photons A/" s (r s ) beyond 
the nonlinear medium can be derived along the follow- 
ing formula considering only photons emitted in photon 



K(t s ) = e m^(u; s )A J dcu t 



(53) 



Using the solutions in Eqs. (|42f and (|43]) we arrive at the 
expression: 

Afsin) = ^ J duj s yfLu~ s J dw' s yftJ s J duji 

Re{exp[i(w s - lj' s )t s ]P*(u s , uj t )F s (uj' s , w<)}. (54) 

Assuming a narrow idler-field spectrum we can rearrange 
the formula in Eq. (|54p as follows: 



JV,(t.) = hio° s J dnReiP^T^n^iT^n)}. (55) 

The functions T s and T' 1 in time domain in Eq. (|55|) 
are derived from their spectral counterparts given in 
Eqs. dSTJ), d45j), and d46j) along the formula: 



P n {T SlTl ) = i- J dtu s J dUiJ^J^^Ui) 

x exp(— iu> s T s ) exp(— iw^Ti), m = s, i (56) 

The number N of coincidence counts caused by a si- 
multaneous detection of both photons from one pair in a 
Hong-Ou-Mandel interferometer is given as: 



x / dti I dt 



'(\r\ A + \t\ A ){(Ei-Hh) 

+ {{rH) 2 {E[r\t 1 )E[+\t 2 )E\-\t 2 -T l ) 
x4 (+) (t 1 -r i ))+c.c.} 



(57) 



r (t) stands for amplitude reflectivity (transmissivity) of 
a beam-splitter in the interferometer and t\ denotes a rel- 
ative time delay introduced between the signal and idler 
photons. Using the solutions in Eqs. (|4"2"| and (I43[) the 
normalized coincidence-count number R n can be derived 
in the form: 



R n {n) = l - p(tz), 



(58) 



where 



x jr s *(cj S) cji)jP(^, Ws )exp[i(w s -uJi)n]j (59) 



and 



Ro = l(\r\ 4 + \t\ 4 ) I dtu s I duj 



xRe 



UJ s UJi 



(60) 
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IV. NONLINEAR LAYERED STRUCTURES 



have to be superposed. So we can write: 



We now consider structures composed of both linear 
and nonlinear layers with thicknesses typically in hun- 
dreds of nm. Fulfilment of phase-matching conditions is 
not important in these short layers, because AkL <C tt 
[L means a typical layer length). Surface SPDC becomes 
important in these structures because volume contribu- 
tions to SPDC from individual layers are weak. A gen- 
eralization of the developed theory to layered structures 
is straightforward and is based on the fact that we de- 
tect only one photon pair. This pair is generated in one 
of the nonlinear layers and propagates to the output of 
the structure. The theory giving volume contributions 
to SPDC has been developed in [23|, H3, HH using a per- 
turbation solution to Schrodinger equation. Here, we re- 
strict ourselves to a scalar model and collinear interaction 
to emphasize the important steps in the derivation. A 
generalization including general directions of fields' prop- 
agation and their polarization properties is straightfor- 
ward following the way presented in 23] . 

First we generalize the momentum operator G; n t writ- 
ten in Eq. ((3|) to include the nonlinear interaction in N 
layers: 



Gi n t(z) 



4e -4 

X/27T 



diOp i du> s I duji 



N 



6{uj p -uj s -uJi)Y^ E V(j 

(=1 a,P,f=F,B 

x 'E^iz^E^^z^E^iz,^) +h. 



(61) 



where the amplitude operators are defined in Zth 

layer: 



2e cAn$(u m ) 



-rect Z! _ 1)Z! (z) 



xa^ } J(z,w m ), m = s,i, a = F,B. (62) 

The pump-field spectral amplitude E^^ (7 = F, B) is 
defined inside the Zth layer similarly as the amplitudes 
in Eq. (|62|) . The nonlinear coefficients ai a p as well as 

indices of refraction (cj m ) characterize the Ith layer 
that extends from zi-\ to zi. Function rect a ,&(.z) equals 
1 for a < z < b and is zero otherwise. 

Photon-pair generation in the Zth layer can be studied 
using the approach and results presented in Sec. II. These 
formulas give us appropriate operators a ma (zi, u> m ) (rn = 
s,i; a — F, B) at the end of the Zth layer. FresnePs 
relations at the boundaries have to be used to 'transfer' 
a generated photon pair outside the boundaries of the 
layered structure. A photon pair can be emitted in any 
of N layers and the corresponding quantum trajectories 



N 

1 = 1 a'=F : B 

m = s, i, a = F, B. (63) 

The coefficients T™f 3 in Eq. §3$ can be derived using 
the propagation matrix method and Fresnel's relations 
at boundaries [13, HI] . 

Properties of photon pairs as described by quantities 
introduced in Sec. Ill and measured by a simultane- 
ous detection of both photons comprising a photon pair 
can be derived from the fourth-order correlation func- 
tion (a^C^ja^^.ja^KOajffwi)). Substituting 
Eq. (|63|) into the definition of the correlation function and 
using the solution for one nonlinear layer as presented in 
Eqs. (gll) and (g3]) we arrive at: 



<a°f , (to' s )aZ K)a^V0<(^i)> = 

/rS.out* / / / \ .7-2, out / \ 



a, a, (3, — F, B. (64) 

The two-photon amplitudes (o — s, i) describing 

a photon pair at the output boundaries with the pho- 
tons propagating in directions indicated by lower indices 
a and (3 can be expressed in terms of two-photon ampli- 
tudes T°J^, characterizing individual layers: 

N 

1=1 a',P'=F,B 

x^^w^Ui), o = s,z; a,(3 = F,B. (65) 

Expressions in Eqs. (143)) and (|4T|) can be rearranged into 
the form: 

f=F,B 

xE(+'Q(w a + Ui) exp[ikg>(u a + lu^Li] 

x exp[-iAk^ a/3 (uj s ,uJi)Li/2} 

xLi smc[Ak^ a/3 (uj s , Ui)L t /2], 

o = s,i; a,/3 — F,B, (66) 



where 



V°'^(w s ,w 2 ) 



,f3 )7 = F,B (67) 



and .9 7 a p{us, w<) = 2idj, a p V^s^t I [2ttc\J (w s )nf (oji)]. 
Superscript (I) denotes quantities appropriate for the 
Ith layer and L[ means the length of this layer 
(Li = zi — Symbol E^'^ occurring in Eq, 
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refers to the spectral amplitude of field p 1 (7 = F, B) 
at the end of the ith layer (at z = z{) and can be 
determined using the propagation matrix formalism. In 



we have assumed that g, 



(0 



a {l) 



deriving Eq. 
and qi l \ P 

The signal-field intensity spectrum S s , its photon flux 
7V S , and coincidence-count interference rate R n in the 
Hong-Ou-Mandel interferometer can then be determined 
using the formulas in Eqs. f52"]l. (T51| , and ([55H5D|) us- 
ing the spectral two-photon amplitudes F s < out and F l < out 
defined in Eq. flBSJ). 

The generalization including fields' propagation under 
nonzero angles of incidence is straightforward following 
the procedure presented above in Sec. II. The key point 
here is that the necessary amplitude corrections assuring 
the fulfilment of continuity requirements for the electric 
and magnetic fields at the boundaries are defined only 
at one side of the boundaries. To be more specific, we 
use nonzero amplitude-operator corrections 6a SF (0,uj s ) 
and Sa SB (0, u> s ) [Sa SF (L, u> s ) and Sa SB (L, u) s )] for the left- 
hand [right-hand] side of the nonlinear medium (layer) 
to describe surface SPDC. When nonzero angles of inci- 
dence are considered, we can argue as follows. The angle 
of incidence of a given field equals (up to the sign) the 
angle of reflection and so their cosines giving multiplica- 
tive factors for projections of electric- and magnetic-field 
amplitudes (lying in the plane of incidence) to the plane 
of boundary are the same. This means that the equa- 
tions written in Eqs. (HHJ), |J3S]), and (33} remain 
valid for any angle of incidence. Equations ([66]) and (|67|) 
can then be used in this case provided that we use the 
z components of wave vectors k p , k s , and ki instead of 
their full lengths. 

Surface SPDC can significantly contribute to the num- 
ber of the generated photon pairs in layered structures. 
As an example, we consider a structure composed of 
25 layers of nonlinear GaN of the thickness of 117 nm 
that sandwich 24 linear layers of A1N of the thickness 
of 180 nm. This structure as a source of photon pairs 
has been studied in detail in [23|, [24[ for s-polarized 
pump (normally-incident), signal, and idler beams and 
the pump wavelength of 664.5 nm. Considering only 
volume contributions efficient photon-pair generation oc- 
curs at degenerate frequencies of the signal and idler 
fields for the signal-field emission angle of 14 deg (see 
Fig. 5 in [23]). The surface contributions lead to ad- 
ditional non-negligible photon-pair generation, as doc- 
umented in Fig. H where the ratio Sf F ' vol+suli /S^ F ' vo1 
of signal-field spectral intensities with ^S FF > vcil + au "^ an( j 
without (S FF,vo1 ) surface contributions is plotted as a 
function of signal- field emission angle 9 S . This ratio is 
minimum under the conditions giving the best construc- 
tive interference (i.e., when ui s =10^/2 and 9 S = 14 deg) 
and equals approximately 2. This means that the sur- 
face contributions roughly double the number of emitted 
photon pairs. It should be emphasized that constructive 
interference between the volume and surface contribu- 
tions plays the key role here. Comparison of the graph 



g FF.vol+surf^g FF.vol 



0)40 

O 

T3 



0.97 0.98 0.99 1.00 1.01 1.02 1.03 
2co /to" 



1 7.48 
10.5 
14.6 
20.5 
28.6 
40.0 



FIG. 4: Contour plot of the ratio S FF ' vol+BUl1 / S FF ' vo1 of the 
signal-field spectra with (Sf i? . vol + surf ) and without (Sf F,vo1 ) 
surface contributions to SPDC as it depends on signal-field 
emission angle 8 3 . Both photons propagate forward. Loga- 
rithmic scale is used on the z axis. 



in Fig. H] with that plotted in Fig. 5 in [23[ indicates that 
the worse the constructive interference inside the layered 
structure the greater the relative contribution of surface 
terms. However, the overall number of generated photon 
pairs is quite low in this region. 

The increase of photon-pair generation rate is caused 
mainly by processes that do not (even roughly) obey 
phase matching conditions. Weights of the surface con- 
tributions for different processes can be judged according 
to the value of parameter V™^g (m = s, i; a, /3, 7 = 
F,B) defined in Eq. ([rJ7| . The following values are 
met in GaN for the studied structure around the point 
where the best constructive interference has been found: 



F,FF ~ V B,BB ~ U - UJ > V F,FB ~ v 



F.BF 



-V 



B.FB 



V 

Vb,bf ~ 2, Vf BB = -Vg >FF » 4. Because the domi- 
nant role in surface SPDC is played by the highly phase- 
mismatched processes, lengths of nonlinear layers have 
to be less or comparable to the coherence length of the 
nonlinear interaction to observe these contributions. The 
ratio N FF ' sur f /JV FF ' vo1 of the numbers of emitted photon 
pairs in a certain spectral region coming from the surface 
(7V FFsurf ) and volume (N FF ' vo1 ) contributions for one 
layer of GaN as plotted in Fig. \5\ indicates that the inclu- 
sion of surface contributions is important for the lengths 
below 1 fim. 

In general, the inclusion of surface contributions leads 
to the broadening of spectral two-photon amplitudes 



J rs - out (w s ,w i ) and T 



z.out i 



,u>i). Thus widths of the 



signal- and idler-field intensity spectra also increase. The 
Schmidt decomposition of two-photon spectral ampli- 
tudes reveals that entropy of its coefficients increases 
provided that the surface terms are included. This 
means that entanglement between the signal and idler 
fields increases due to surface SPDC. On the other hand 
and considering pulsed pumping, photon fluxes M and 
coincidence-count interference rates R n in a Hong-Ou- 
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L (I0" 6 m) 

FIG. 5: Ratio N FF ' mTt /N FF ' vo1 of the numbers of emit- 
ted photon pairs originating in surface (jV FF ' surf ) and vol- 
ume (jV FF,vo1 ) down-conversion from one layer of GaN of 
length L pumped at X p = 664.5 nm at normal incidence. 
Pairs with frequencies in the interval Ail used in the graph 
in Fig. U and propagating along the +z axis are collected; 
N — J An diOs J AQ duii n(uj s , OJi) and n is given in Eq. (|49}. 
Logarithmic scale on the y axis is used. 

0.15 -i 





FIG. 7: Normalized coincidence-count rate R„ in Hong-Ou- 
Mandel interferometer as a function of relative time delay r; 
is shown provided that volume (solid line with *), surface 
(solid line with A), and volume + surface (solid line with 
o) contributions are included. The structure is pumped by a 
Gaussian pulse 250 fs long. 



continuities is present inside the nonlinear crystal. The 
surface effects are also expected in wave-guiding geome- 
tries where the conditions of total reflection of fields at 
the boundaries are met. This might be important mainly 
in photonic fibers. The studied surface effects are by no 
means restricted to ID geometries, even stronger surface 
contributions might arise from 2D or 3D nonlinear struc- 
tures. An effective enhancement of the nonlinear interac- 
tion caused by the surface effects should also be observed 
in stimulated processes like second-harmonic gener- 
ation. 



FIG. 6: Signal-field photon flux A/" FF including volume (solid 
line with *), surface (solid line with A), and volume + surface 
(solid line with o) contributions. The structure is pumped by 
a Gaussian pulse with the duration of 250 fs (for details, see 

Mandel interferometer as temporal characteristics be- 
come narrower as documented in Figs. [6]and[7] The over- 
all photon flux Af of, e.g., the pulsed signal field occurs 
earlier at the output of the nonlinear structure compared 
to the case of only the volume interaction because the 
structure is 'less- resonant' ('less-transparent') for the sur- 
face contributions than for the volume ones. Phase mod- 
ulation of the two-photon spectral amplitude character- 
izing the surface contributions provides the coincidence- 
count pattern R n (ri) in the form of a global dip with 
visibility equal to 1 and two small side-dips (see Fig. [7]). 

Surface effects in SPDC occur at discontinuities of 
X^ 2 ) nonlinearity. Such discontinuities are found also 
in periodically-poled nonlinear crystals like LiNbC>3. In 
these crystals, an enhancement of photon-pair generation 
rates is expected provided that the period of poling is suf- 
ficiently short so that a sufficiency large number of dis- 



V. CONCLUSIONS 

Surface spontaneous parametric down-conversion has 
been predicted combining the solution of quantum 
Heisenberg equations and the continuity requirements at 
boundaries. Formulas for the determination of the num- 
ber of generated photon pairs, spectra, and photon-fluxes 
of the down-converted fields as well as coincidence-count 
rates in a Hong-Ou-Mandel interferometer have been de- 
rived using generalized signal- and idler-field two-photon 
spectral amplitudes. It has been shown that surface con- 
tributions from the input and output boundaries of a 
nonlinear crystal give structures similar to those charac- 
terizing the volume contributions. Nevertheless they are 
weak. The surface contributions are important when- 
ever strongly phase-mismatched nonlinear interactions 
give considerable contributions, e.g., in nonlinear layered 
structures. An example of a GaN/AIN structure com- 
posed of several tens of layers has shown that the surface 
and volume contributions can be comparable in their am- 
plitudes. This shows that the role of surface effects in 
other nonlinear structures like periodically-poled materi- 
als, nonlinear wave-guiding structures or structures with 
stimulated processes should be addressed. 
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